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Abstract 

A problem concerning the shift of roots of a system of algebraic equations is investigated. Its conser- 
vation and decomposition of a multiple root into simple roots are discussed. 

1 Introduction 

The central subject of the present paper is an investigation on the shift of roots of the system of algebraic 
equations. Our central Theorem 13 . II states that the number of real roots of a system located in a compact 
set does not change after a sufficiently small perturbation of the system. 

As a matter of fact, this kind of fact has been well known to those who study the deformation of the 
singularities of diffcrcntiable mappings. It is, however, a non trivial question how small this perturbation 
shall be so that the number of roots in a given compact set remains unchanged. All the ever existing theorems 
(see §12.6) do not precise the size of the compact set and the perturbation of the system under question. 
They state simply that for a compact set and perturbation, both of them small enough, the invariance of 
the number of roots holds. This situation can be explained by the fact that they simply treat the notion of 
local algebra, and consequently they are valid only in the germ sense. We try to give here an estimate on 
the size of admissible perturbation of the system for a fixed compact set. 

Furthermore, we give a result about the decomposition of multiple roots into simple roots. In particular, 
our Theorem 14.21 assures us the existence of a deformed system of the original system that possesses only 
simple roots. 

The paper is organized as follows. In §2 preliminary results are presented. In §3 we state that a slightly 
deformed system has the same number of zeros as the original system in taking the multiplicities into account. 
In §4 we give a result about the decomposition of multiple roots into simple roots. The paper ends in §5 
with some concluding remarks. 



2 Preliminary results 



Let us consider the following system of algebraic equations with real coefficients Oq 

fi(xi,x 2 ,...,x n ) = 0, 
h{xx,x 2 , ...,x n ) = 0, 



fn{xi,X 2 , ...,X n )=0, 



w. 



(1) 
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with 



fi(x\, X2, ■ ■ ■ , x n ) 



J l ^2 

loci = oil + a 2 + ■ ■ ■ + a n , 



\a\—Trii 



(2) 



where the degrees of polynomials are ordered as follows: 

771 1 1712 ^ • • ■ m n . 

Let us consider the situation where the gradient ideal (^$^fe, ■ ■ ■ , af - /^ over lRL[x] contains certain power 
of maximal ideal m k . That is to say, 



dfi(x) dfi(x) dfiix) 



dx\ ' 8x2 



dx r , 



_j \a 1 x 2 x n / QiH \. an ^k ■ 



(3) 



Let us note by M a (x) — x" 1 ■ ■ ■ x" n , a\ + • • • + a n = k. It is possible to consider the set of such 
monomials as a basis of m fc . The dimension /i„(fc) of the basis of the ideal m fc can be calculated by the 
following recurrent relation: 





k 




i=0 




= k+l, 


A*3(fc) 


_ (k + 2 



Evidently, /i ra (fc) is the number of the entire lattice points on an (n — 1) -dimensional face of the n-simplex: 

= # {("i, ■ • ■ ,a n ) € Z| : a\ H h a„ = fc} . 

We recall the notation that: 

supp (f(x) = {a e Z n ; <yS Q ^0}, 

for a Laurent polynomial tp(x) = L PaX a . We call a germ y?(x) convenient at zero when the Newton 

diagram of it at zero contains non-compact part of all coordinate axes (cf. |2])- In other words (p(x) is 
convenient at zero if it admits the representation, 

n 
i=l 

for ^ ^ 1 and a certain polynomial R(x). It is easy to see that if fe(x) has a convenient germ at zero then 
their exists fc^l such that the condition is satisfied. 
Suppose that a polynomial vector: 



/Mx)\ 

¥>2(X) 

\ip n {x)J 



= F ■ 







V J 
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where F is an invertible constant matrix and f{x) £ m fc+1 and degip(x) = k' ^ k + 1, i.e. supp f{x) C {a £ 
Z™ ; k + 1 ^ |a| ^ k'}. The question we pose concerns the behavior of roots of a system: 



f 1 (x)+t<p 1 (x)=0, 
f 2 (x)+tcp 2 (x)=0, 

f n (x) + tip n (x) = 0, 



with t £ [0, 1] C t as a parameter. 

To formulate further statements in a proper way, we introduce notations: 



(4) 



where the vectors v 



fi(x) = afx^ + afx^ +■■■ + og: 
( u j*i> v j% ••■■> v j%) > 1 < 3 < A *> satisfy: 



(5) 



((1, • ■ • , 1), *f ) = ((1, • ■ ■ , 1), 4 l) ) = ■ ■ ■ = ((1, • • • , 1), - m,. 

In general, it is not easy to formulate a sufficient condition on /e(x) so that the condition |J3J| holds. We 
propose here a simple necessary condition for that. 

Propositin 2.1 The following isomorphism |3J) is necessary so that the condition |3J) holds, 

lattice spanned over Z of ju^ — t^, ■ • ■ , vf^ — — 

{5eZ";((l,.,l),a) =0}. 

Proof If the condition docs not hold, it is evidently impossible to create all monomials x a with 
|a| = k as a linear combination of ■£^-fi(x)'s. Q.E.D. 

Remark 2.1 This isomorphism can be realized by shifting each lattice point v of the right hand side of 
Relation |5|) towards the lattice point v + (m„ 0, . . . , 0). 

We give now an example for which condition does not hold. 

Example 2.1 Let us consider the following system: 

fx(xi,x 2 ) = xf + a^xfxl + a-px% = 0, 



h{xx,x 2 ) = a ( ^x\ 2 + a { 2 2) x\xt + a% ] x\ 2 = 0, 



(6) 



whe 



4<4 <) 4 i) ^0, i = 1,2. 



For these polynomials the lattice define on the left hand side of Relation J3|) is isomorphic to 



a£Z 2 ; ((l,l),a)=0; ^£Z 2 



Thus, in this case, the condition |3J) does not hold. 



From now on, we use the notation fe(x, a) instead of fe(x) if we want to emphasize its dependence on the 
coefficients a = (a\ , ■ ■ ■ , a l x ). For a set of polynomials A\{x, a), A 2 (x, a), . . . , A T V (x, a) £ R[x, a] homogeneous 
in variables a we consider the following linear combinations: 

A^ +1 (x, a)=J2 j( j \a)xW A\{x, a), j = 1,2, . . . K+1 , 



where £ Z> and t| (a) are linear polynomials in variables a. 
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Propositin 2.2 Let as consider the chain of polynomials sets: 

{a[(x, a),A r 2 (x, a),..., A r „ r (x, a)} , r = 0, 1, 2, . . . 

as above with 

d d 

Suppose that for certain r — L, some of A%(x,a) 's coincides with X^(a)M s (x). That is to say there exists 
h[ s \x,a), . . . ,hn\x, a) e R[x,a] such that 

" fl 
X^(a)M s (x) = a fe^' 

T/ien deg a (a) = L and deg a (x, a) = L - 1 . 

Proof After the definition of the recursive process to create A^ +1 (x,a) from A r (x,a) it is clear that 
Aj(x, a) is a homogeneous polynomial of degree r in a. The statement is the direct consequence of this fact. 
Q.E.D. 

Example 2.2 We consider the following example 

fe{x 1 ,x 2 ,a) = f(x,a) = a^x\ +a 2 xjxl + a 3 xl = 0, 

with 

d 

- — f(x, a) = 5aixj + 2a 2 xix\, 
ox i 

d 

— — fix, a) — "ia 2 x\x\ + 5a 3 x 2 . 
ox 2 

Then we have the following chain of polynomials to getx a2 — x\x\ as a linear combination of0\ = -^f{x,a) 

d 

8X2 



and 2 = J-f{x,a): 



A\(x,a) 
A\{x,a) 
A\{x,a) 



= (2a 2 x\x 2 — haix x x 2 ) 6i, 
= (3a 2 xf - 5a3xfx%) 9 2 , 
= 5a 3 xlx 2 2 , 



Thus we have: 



A\(x, a) := 3a 2 A\(x, a), 

A 2 (x,a) := 3a 2 A 2 (x,a) + 5a 3 Al(x,a), 

Af(x,a) := 3a 2 Al(x,a), 
A\{x, a) := 5aiA 2 (x, a), 

A\{x, a) := 3a 2 Al(x, a) + 5aiA^(x, a) = (2 2 3 3 a^ + 5 5 a 2 a%)x\ x\. 



h^\x,a) = (3a2) 3 (2a2Xia;2 — daix^x?,), 

h 2 2 \x,a) — (5ai) 2 ((3a 2 ) 2 a;^ — lha 2 azx\x 2 — (5a 3 ) 2 xfx 2 ). 
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For the case of x a5 = x\x\, we have: 

\( 5 )(a) = 5oi(2 2 3 3 a^ + 5 5 a 2 a|), 

hi(x, a) — ba^a^x\x\ + 2 3 3 Baia^xfx 2 , 

h^\x,a) = —2 • b 2 a\a2((3a2) 2 x\ — \ha2a?,x\x\ ~ (haz) 2 x\x^)- 

For the case of x ai ° = x 2 ° j we have: 

X^°\a) = 5a 3 (2 2 3 3 al + 5 5 a 2 a 3 2 ), 
h^ 0) (x,a) = (2 2 3 3 a| + 5 5 a 2 ai)zf - 

— 3a2(5ai) 2 ((3a2) 2 x^ — lba2a^x\x\ — (5as) x\x^) , 
h 2 W ^ (x, a) — — (3a2) 4 (2a2Xixl — baix^x^)- 



3 The number of roots of a deformed system 

In this section we state that a slightly deformed system has the same number of zeros as the original system 
in taking the multiplicities into account. We recall here that the index ^ £ [1, n] has been fixed so that fi{x) 
satisfies the condition 10. 



Definition 3.1 We introduce the norm: 

IMI = E l«l 

a£ supp ipi (x) 

where <pg(x) = J2 a <Esu PPVe (x) L Pn,aX a ■ We name the following value by C(a): 

C(a) = max | max h{ s }(x,a)x^ 

\1^<" |/3|<fc'-fc-l 

for some compact set K and hj S J(x,a) — jj7rn~r after the notation of Proposition 

Remark 3.1 In general we can not give any reasonable estimate on C(a). In the above Example 
hjj(x,a) contains coefficients of the form: 

polynomial of degree 5 in (a\, 02, 03) 
a 3 (2 2 3 3 a| + 5 5 a 2 a 3 ) ' 

This value can be as large as possible if the denominator is very near to zero. The coefficients of h 3 s i(x,a) 
contain rational functions in the variable a, with denominators \( s '(a) introduced in Provosition \2/A □ 

Before formulating our main theorem, we recall a simple lemma of linear algebra. 

Lemma 3.1 Let us consider /1 x \x matrix A — (a.y) £ End(R M ). If \a%j\ < -pr, then (id M + A) is invertible. 

Proof By straightforward calculation of the determinant of (id M + A) we have: 

det(id M +A) = l + a n +a 22 -\ h + R(a), 

where R(a) is a polynomial containing (/i 2 — /x — 1) terms of monomials in (cty) whose degrees are higher 
than or equal to two and less than or equal to /1. Evidently, under the condition \a%j\ < A? we obtain that 
det(id M + A) ^ 0. Q.E.D. 
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Theorem 3.1 The number of simple roots of the system ^ inside of a compact set K coincides with that 
of the system if t satisfies the following inequality: 



t < 



\W\\C{a)tf 



(7) 



Proof Our strategy consists in the construction of a homotopy that connects the simple roots of system Q 
and those of 

Suppose that we succeed in constructing a homotopy x(t), ^ r ^ t with x(0) = x such that 

/ s (x(r)) + Tip s (x(r)) = f s (x), 1 < s < n, 
then the vector field along it satisfies the following equality: 

/A(z(t))\ 

71 / f-^f^W 



d_ 



1 Mr)) » 



V/n(x(r))/ 

In applying this relation to system Q, we get. 

r //iWt))\ 

I \f n (x(r))J 
Further, we shall realize a smooth homotopy 



tF 



/¥>(*) \ 





V o J 



f2(x(r)) 
V/„(x(r))7 



F 



/<p(x)\ /0\ 
1 

V o ) 



We remember that we denoted the basis of m k by M ai (x), 1 ^ i ^ /i = /i„(fc). The condition 10 entails the 
following relation: 

/M ai (x) \ 



M Q2 (x) 



W„(i)/ 



a/2 

5a: i 


+ i j(2) 


/Mx 

a/2 

dx 2 









+ • • • + 



a/2 1 
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for some polynomial entry rank-1 (fj, x n) matrices H' 1 ' , H^ 2 ', . . . , of the form: 



(8) 



r /o 



••• o\ 



••• ••• 



I Vo ••■ htt o ... oy 



h[ s) (x,a) 



where h^j(x,a) — "j^)^" -1 after the notation of Proposition 12.21 concentrated at the ^-th column of the 
matrix ifW. One rewrites the relation JSJ) as follows: 



/M ai (x)\ 
M a2 (x) 

\M a (x)J 







J 2 + T<^ 2 
V /n + T(p n / 
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M*)\ 




v o y 



(9) 



As we supposed that (p(x) € m fc+1 , it is easy to see that: 







= A{x) 



V / 

with certain polynomial (fi x /i) matrix A(x), where: 

a(x) = ; 

By recalling ©, we obtain an equation as follows: 



(M ai {x)\ 
M a2 (x) 



(10) 



+ tA(i)) 



(M ai (x)\ 
M a2 (x) » 



9a;,; 



Supposing that r is very small, get the inverse to: 

(id M + Tyl(x)), 



J2 + T(f 2 
V /« + Tip n ) 



(11) 



in the domain {x; det(irf A1 + tA(i)) ^ 0}. The inequality J7J ensures the invertibility of the matrix To 
show this, in view of Lemma 13. II it is enough to verify that for such a value of r we have: 



r | max max 



< 



In other words, it is enough to prove that: 



max max 



<|M| C(a). 



(12) 



(13) 



We remember that suppgi^(x) c{a£ Z n ; fc' — m n < \a\ ^ fe' — mj. This is a direct consequence of lfTU)l . 
As we have suppipe(x) c {a € Z"; fc + 1 ^ |a| < fc'} , we can find for every l^A^na series of polynomials 
#2 (*),..., (sc) such that: 



In terms of these polynomials: 



A w ( x ) = ^ 5 g(,) Mc (4 



9x A 



(14) 



and 



supp^l(x) c{ae Z"; ^ \a\ < jfe' - k - 1} , 



7 



the absolute value of each coefficient of Q\{x) can be estimated by \\ip\\ after (|14|l above. We replace Q c ^(x) 



by IMI x (e^, 



>--*»i>pe&w •'' and we get the ;n( ' l| '' ; ' :i,v: 



max max 



g<P{x) < [ J2 Hl^,a| 

aE supp tp(x) 



max N 

i«j<n ^ eUr=iSupp ^£M VU M =iSuppMs(:l;) 



E 



max 



(1,0)0^ 



where A \ B — {a — (3 G Z> ;a! £ A, [3 E B}. The Relation l|14|) explains the summand of the above 
inequality. Therefore, if we set C{a) as in Definition 13.11 we obtain the inequality (|13|l . Evidently, C(a) 
depends not on the coefficients of ipi(x), . . . , ip n (x) but on the powers k and k' . This proves the invertibility 
of the matrix iJUJ. Thus, 



/M ai {x)\ 
M a2 (x) 



= (^ + rA(,)r 1 (^i/« — 



\M a (x)J 



I fl + Tifl \ 
f-2 + Tlfi 2 

V /„ + Tip n I 



On the other hand: 



tP2(x) 



F ■ 



\ip n (x)J 

for some rank-1 (n x /i) polynomial matrix: 



/<p(x)\ 




V J 



F-G 



/M ai (x)\ 
M a2 (x) 

\M a (x)J 



G = 



( 9i(x) g2(x) 




\ 

If we apply F ■ G from the left to the relation (|1 5|) . we get: 







/ 



\ip n {x)J 



= F-G{id^+TA{x)y 1 [Y J H {i) ^ 



or 







V / 

This relation gives rise to an inequality: 



^ n 





v J 





Vo 



Vi(x,r)Q 









0/ 



_d_ 

dxi 



( fl + Tip! \ 
\fn+ TLp n I 

V /„ + T^ n / 
/ A + T(fl \ 



/ 2 + Tip 2 
V /n + T^n 



(15) 



(16) 



(17) 



(18) 
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where the matrix in front of the derivative has single non zero £-th column. That is to say we obtain the 
following equalities: 



<p(x) 



Vi (x, t) ~a~~ ) (ft+T(pt), 



dx 



vi=i ^ / (19) 

ipi(x) = F ei ip(x) 

The estimate 10 ensures that Vi(x,r) are real analytic in x € K. Thus we have constructed a vector field 
corresponding to the homotopy we need. Q.E.D. 

Example 3.1 Let us consider the following system ^: 

fx(xi,x 2 ) = x\ - x\ - 1 = 0, 



f 2 {xi,x 2 ) = x 2 +xj-2 = 0. 

This system has four real solutions within the square [— 2,2] 2 : 

(±VT5,±\/(L5) » (±1.22474487139159, ±0.70710678118655). 
// we perturb this system with a cubic monomial <p(x) = X\x\ as follows: 

Fi(xi, x 2 ) = x\ — x\ + tx\x\ — 1=0, 
F 2 (x!,x 2 ) = x\ +x 2 2 - 2 = 0, 



(20) 



(21) 



we calculate the constants C(a) — |, \\ip\\ = 3,/i = 2. Therefore we have four solutions of the system l^21\) if 
t < j^E' ^ n P ar ^ cu ^ ar ! if we use the value t — 0.033 < 1/30 by applying the rootfinding method of llttf we 
obtain the following four solutions: 

( 1.22054232589618, ±0.71433635683474), 
(-1.22879457180552, ±0.70004564158438). 

Furthermore, if we perturb both equations of the system \2ty with the same cubic monomial tp(x) = x\x\ as 
follows: 

Fi(xx,x 2 ) = x\ - x\ + tx\x\ -1=0, 

(22) 

F 2 (xi,x 2 ) = x\ + x\ + txixl -2 = 0, 

we again calculate the constants C(a) — I, \\<p\\ = 3,/i = 2. Therefore we have four solutions of the 
system l^22\) if t < 2 * 5 . In particular, if we use the value t — 0.033 < 1/30 by applying the rootfinding 
method of we obtain the following four solutions: 

( 1.21652265747566, ±0.70710678118655), 
(-1.23302265747566, ±0.70710678118655). 



Example 3.2 Let us consider the following system: 



fi(x 1 ,x 2 ,x 3 ) = \x\ + ^ 



81 



0, 



f 2 {xi,x 2 ,X3) = [xx + x 2 ) 2 + 36(xi - x 2 ) 2 - 9x1 = 0, 



2 2 

X-t n Xr> 



f 3 (x u x 2 ,x 3 ) = ^+xl + -^--l = 0. 



(23) 
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By applying the rootfinding method of JflSf we obtain the following sixteen real solutions within the cube 
[-2,2] 3 : 

0.62830967308983, 0.91412675198426, ±0.76883755100759), 
-0.62830967308983, -0.91412675198426, ±0.76883755100759), 
0.49635596537865, 0.91441703848857, ±0.95929271740718), 
-0.49635596537865, -0.91441703848857, ±0.95929271740718), 
1.11731818404380, 0.76796989195429, ±0.93973420474984), 
-1.11731818404380, -0.76796989195429, ±0.93973420474984), 
1.22450432822695, 0.66467487192937, ±1.28459776563576), 
(-1.22450432822695, -0.66467487192937, ±1.28459776563576). 

We observe that these roots are invariant under the actions of a group G := (Z/2Z) x (Z/2Z) generated 
by two generators (x\, X2, x 3 ) > (— Xi, — X2, Xs) and (x±,X2,x 3 ) <— > (x\, X2, — £3) due to the invariance of 
the system $2,1\) itself under the same group action. If we perturb this system with a quadratic monomial 
if(x) = x\ as follows: 



fi(xi,x 2 ,x 3 ) = [x( + — 



f2(xi,x 2 ,x 3 ) = (xi + x 2 ) 2 + 36{xi - x 2 ) 2 



9x 2 3 



f 3 {x 1 ,x 2 ,x 3 ) 



1 = 0, 



--^=0 
81 ' 

tx\ = 0, 



(24) 



we calculate the constants C(a) — i, \\(p\\ = 2, /i = 3. Therefore we have sixteen solutions of the system 

if t < ^j-. In particular, if we use the value t = 0.1 < 1/9 by applying the rootfinding method of ]13f we obtain 

the following sixteen real solutions within the cube [— 2,2] 3 : 

0.63087661393950, 0.91351892559324, ±0.77060795720733), 
-0.63087661393950, -0.91351892559324, ±0.77060795720733), 
0.49896002229193, 0.91405620623649, ±0.95934810309529), 
-0.49896002229193, -0.91405620623649, ±0.95934810309529), 
1.11568183127565, 0.76857206607484, ±0.93967783245553), 
-1.11568183127565, -0.76857206607484, ±0.93967783245553), 
1.22357424633595, 0.66527556809517, ±1.28379297777855), 
-1.22357424633595, -0.66527556809517, ±1.28379297777855). 

One remarks here also the invariance of the roots under the above mentioned group G due to the invariance 
of the system ji^| ) itself. □ 

As for the equation we establish the following theorem: 

Theorem 3.2 Let us consider a system of algebraic equations obtained as a perturbation of f7|) : 

Fi(xi,X2, ■ ■ ■ ,x n ) = 0, 
F 2 {x ll x 2l ■■■,x n ) 



Suppose that on a ball B r = {i£ 



0. 



F n (x 1 ,x 2 ,...,x n ) = 0. 
x\ ^ r} we have: 



(25) 



rank I -^—fk{x 
1 oxj 



(26) 
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Furthermore, we impose a condition on (Fi, F2, ■ ■ . , F n ): 



max 



vector component of 



CPi-/i)(z) 
(F2-f 2 )(x) 



max 



vector component of ( — — -F k {y) 





~ fn)(x 


(Fi 




(F 2 





< 



(27) 



(F n -f n )(x) 



< E, 



(28) 



where e is strictly less than distance of any root of in B r io the boundary 9B r . Suppose that the 
system 0) /ias no multiple real roots. 
Under these assumption the equality: 



#{real simple roots of Q) inB r } = #{reoi simple roots of mB r } 



holds. 



Proof We solve the homotopy equation with respect to smooth diffeomorphism 

Xi + hi(x, t), ^ i ^ n, 

that satisfies 

hi(x,0) = 0, 

f k (x + h (x,T), ...,x n + h n (x,r)) = f k (x , ...,£„) + 

+T(F k (x) - f k (x)), 0<r<l. 

The system gives rise to a system of {n + 1) nonlinear differential equations: 

9ft 



^2 -Q^-djf k (x + ho(x,r), ...,x n + h n (x,r)) = F k {x) - f k (x), 



3=0 



(29) 
(30) 

(31) 



for k = 0, 1, . . . , n. From the assumption (|26f) . Eq. I|31|l is always solvable in the class of real analytic functions 
so far as 



^0 



after Cauchy-Kovalevskaya's theorem on the quasidinear partial differential equation. After the condi- 
tions (|2"?jl - l(2"%|l and Eq. IpTTJ) . is always strictly less than e. Therefore |/ij(r)| < er and |i + hj(r)\ < 
\x\ + e t < r for x root of located in the ball B r . Thus the homotopy equation admits a real analytic 
solution that connects x £ B r with x + h(x, 1) € B r . Q.E.D. 

Corollary 3.1 If the conditions on the analyticity of the homotopy constructed in the above Theorems \'j.l\ 
and \3.<H are fulfilled then none of the roots of the deformed system f35p crashes with another and consequently 
no new multiple roots are created after the proposed perturbation. 

Proof Assume that after the proposed perturbation a multiple root is created. Then the homotopy 
constructed in the above Theorems 13.11 and 13.21 looses its analyticity with respect to the parameter r. 
Q.E.D. 
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4 Decomposition of multiple roots 

In this section we recall facts about the decomposition of multiple roots into simple roots. 
Definition 4.1 For the system fl)) . we define the Jacobian function: 



(j/)(x)=det 



dxi dx n 



dfn . _ Qfn 

dx\ dx n 



(32) 



Let us denote by Vjf — {x £ R™ : jf(x) = 0} the zero set of jf(x). We use the notation Qjf for the following 
set 

V 3f n{x:{fl + --- + f 2 n ){x) = {)}. 

Then we have the following result: 

Theorem 4.1 If Qjf = then the system £7J) has only simple roots, while if Qjf ^ then the system 0J) 
has multiple roots. 

Proof The proof follows from well known facts in the singularity theory [5] . We prove the contrapositive 
of the statements. Namely we can easily see that the existence of multiple roots yields non-emptiness of 
Qjf. On the other hand, the existence of only simple roots entails the emptiness of it Qjf. Q.E.D. 

Definition 4.2 Let us denote by (J/) a vector valued ideal 

dh ... dfi_ 

dxi dx n 



(Jf) 



(33) 



dfn . df n 

. dxi dx n 

Theorem 4.2 Let us consider a system like for which we know that it possesses m real roots with 
multiplicities n%,n2, ■■ ■ ,n m where 1 ^ nj for 1 ^ j ^ m. Then there exists a vector polynomial: 



e (R[x]) n / (Jf) 



such that the system of equations 



H n {x) 

(f 1 +H 1 )(x)=0, 
(/„ + #„) Or) = 0, 



(34) 



(35) 



has n\ + 71,2 



n m simple real roots. 



Proof Wc remark that 



(f l + H l ){x)=Y J Kk(x')x? 



k=t) 



with hi^ mi (x') ^ Q,hifi(x') ^ 0, x' — (x2, • • • , x n ), after certain permutation of variables x. It is well known 
that, there exists a perturbation Hi(x) such that the equation X)fc=o ^ l i.k{x')x" li ^ k — has m, simple roots 
for a codimension 1 set of x' (cf. |2])-This fact entails that the system (|35|l also possesses as much simple 
roots as Q has. Q.E.D. 
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Remark 4.1 One can understand this theorem by an intuitive way. Let us denote by Ii :— supp(fi + Hi) 
the set of powers present in the polynomial fi + Hi. If the discriminant of the system 



(fi+Hi)(x)= ^/ a /, 



(36) 



aeli 



i ij.J does not vanish, then the roots of the system IHb]) are all simple. That is 



say, A(/ Q ) S 

to say the set of the coefficients of the system if#6)) for which the system has multiple roots is of codimension 
one in the space of coefficients R^>=i This fact is known under the name of Bertini-Sard theorem fity . 



Remark 4.2 It is worthy to notice that one shall choose a proper vector polynomial J^J), to get distinct 
simple roots for the deformed system $35}) . For example if f\ — x 3 the fi — has a triple root at x — 0. // 



we take H\ 



-x then fi 

x 



Ht 



x = x(x — 1) = has 3 distinct roots x = — 1, 0, 1, while for H\ 



the equation x 3 + x = has 1 real simple root at x — and two distinct complex roots. 
Example 4.1 Let us consider the following system: 

fi(xi,x 2 ) =xl-x%-l = 0, 



f2(xi,x 2 ) = xf + : 



1 = 0. 



(37) 



This system has two multiple real solutions (±1,0) within the square [— 2,2] 2 . 

// we perturb this system with H\ = and the simple linear polynomial H 2 — t{x\ — 2) where < t ^ 0.5 
as follows: 



(fi + H l ){x l ,x 2 ) 



x\ 



1 = 0, 



{H + H 2 ){x 1 , x 2 ) = x\ + x 2 2 + t{x x - 2) - 1 = 0, 



(38) 



then we have four simple real solutions. In particular, if we use the value t = 0.5 by applying the rootfinding 
method of 1131/ we obtain the following four real simple solutions: 

(1.07123233675477, ±0.38410769233261), 

(-1.20970135357686, ±0.68071827127359). 
While if we use the value t — 0.025 we obtain the following four real simple solutions: 

(1.00412951827050, ±0.09097301502177), 

(-1.01237171332486, ±0.15778620326351). 
Finally, if we use the value t = 0.0125 we obtain the following four real simple solutions: 

(1.00207398824224, ±0.06443817123186), 

(- 1 .00621769007449, ±0. 11168724107454) . 

Example 4.2 Let us consider the following system: 



fl(X!,X 2 ,X 3 ) 



f 2 (xi,x 2 ,x 3 ) 
f 3 (x 1 ,x 2 ,x 3 ) 



((xi - x 2 ) 3 - x 2 3 (xi +x 2 ~ x 3 )) 
((-Xi - x 2 ) 3 - x 2 3 {-xi +x 2 - x 3 )) 
((-xi + x 2 ) 3 - xl(-xi - x 2 - x 3 )) 
((xi + x 2 ) 3 ~ xl(xi -x 2 - x 3 )) = 0, 



2 2 



+ -Z - 1 = 0. 



(39) 
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This system has eight simple real solutions: 

(±0.25926718242254, ±1.21300057180546, ±1.75418919109753), 
and eight triple real solutions: 

(±0.68824720161168, ±0.68824720161168, ±1.37649440322337). 

within the cube [—2, 2] 3 . We observe that these roots are invariant under the actions of a group T := (Z/2Z) x 
(Z/2Z) x (Z/2Z) generated by three generators (xi, X2,X3) i— ► (— xi, X2, X3), (xi, X2, X3) 1— ► (x\, —X2, X3) and 
(xi, X2, X3) 1— ► (xi, X2, — X3) due to the invariance of the system 1 ,9.9)) itself under the same group action. 
If we perturb this system with H\ where < t ^ 0.5 and H2 = H3 = as follows: 

{fi + H 1 )(x 1 ,x 2 ,x 3 ) = 

((xi - x 2 ) 3 - t{x\ - x 2 ) - xl(xi +X2 - X3)) 

((-xx - x 2 ) 3 - t(-xi - x 2 ) - x§(-a;i +x 2 - x 3 )) 

((-£Ci + x 2 ) 3 - t(-Xi + x 2 ) - xl(-xi - x 2 - x 3 )) 

((xx + x 2 ) 3 - t( Xl + x 2 ) - xl( Xl - x 2 - x 3 )) =0, (40) 

x 2 

{f 2 +H 2 )(xx,X2,X3) = xl+xl--2- = 0, 

2 2 

(f 3 +H 3 )(x 1 ,x 2l x 3 ) =x 2 x + ^ + ^--l = 0, 

then we have thirty two simple real solutions. In particular, if we use the value t — 0.5 we obtain the following 
eight real solutions which are shifts of the simple solutions to System \3ytl : 

(±0.27142016486929, ±1.20645760731621, ±1.74883324771051). 
Also we obtain the following twenty four simple real solution: 

(±0.68824720161168, ±0.68824720161168, ±1.37649440322337), 
(±0.78897550317143, ±0.32932116069209, ±1.20907797224513), 
(±0.44474589932680, ±1.07278013064881, ±1.64234961179579), 

that originate from the triple solutions to System ftfyp . We remark that the first ones of the above solutions 
coincide with the triple solutions to the original system. These roots are also invariant under actions of the 
group r. □ 

5 Concluding remarks 

A problem concerning the shift of roots of a system of algebraic equations has been investigated. Its 
conservation and decomposition of a multiple root into simple roots have been discussed. 

To this end, with our central Theorem 13.11 we show that the number of real roots of a system located 
in a compact set does not change after a sufficiently small perturbation of the system. This theorem 
can be applied to high dimensional CAD where it is sometimes needed to calculate intersection of several 
hypersurfaces that are perturbation of a set of original (unperturbed) hypersurfaces. For example, to draw 
a 3D (three dimensional) picture of a real algebraic surface that obtained as a deformation of a known one, 
the question of the shift of roots plays essential role §9.6]. We hope that our results to this direction 
represent certain interests to those who are interested in the application of algebraic equations to computer 
graphics. 

Furthermore, we give a result about the decomposition of multiple roots into simple roots. In particular, 
our Theorem 14 . 21 assures us the existence of a deformed system ff 3511 of the original system QJ that possesses 
only simple roots. This result can be used in many cases including the computation of the topological 
degree ^ElEl^nEl i n order to examine the solution set of a system of equations and to obtain information 
on the existence of solutions, their number and their nature ^El 101 HI HH • 
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